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SPANNING FORESTS AND THE VECTOR BUNDLE LAPLACIAN^ 

By Richard Kenyon 
Brown University 

Dedicated to the memory of Oded Schramm 

The classical matrix-tree theorem relates the determinant of the 
combinatorial Laplacian on a graph to the number of spanning trees. 
We generalize this result to Laplacians on one- and two-dimensional 
vector bundles, giving a combinatorial interpretation of their deter- 
minants in terms of so-called cycle rooted spanning forests (CRSFs). 
We construct natural measures on CRSFs for which the edges form 
a determinantal process. 

This theory gives a natural generalization of the spanning tree 
process adapted to graphs embedded on surfaces. We give a number 
of other applications, for example, we compute the probability that 
a loop-erased random walk on a planar graph between two vertices 
on the outer boundary passes left of two given faces. This probability 
cannot be computed using the standard Laplacian alone. 

1. Introduction. The classical matrix-tree theorem, which is usually at- 
tributed to either Kirchhoff [15] or Brooks et al. [2], states that the product of 
the nonzero eigenvalues of the combinatorial Laplacian on a connected graph 
is equal to the number of rooted spanning trees of that graph. This theorem 
has extensions to graphs with weighted edges (resistor networks) and more 
generally to Markov chains. Forman [5] extended the theorem in a more in- 
teresting direction, to the setting of a line bundle on a graph. In this paper, 
we reprove his result and extend the theorem further to two-dimensional 
vector bundles on graphs (with SL2C connection). Our main theorem gives 
a combinatorial interpretation of the determinant of the Laplacian as a sum 
over cycle-rooted spanning forests (CRSFs). These are simply collections of 
edges each of whose connected components has as many vertices as edges 



Received January 2010; revised April 2010. 

^Supported by the NSF. 

AMS 2000 subject classifications. 82B20. 

Key words and phrases. Discrete Laplacian, quaternion, spanning tree, resistor net- 
work. 

This is an electronic reprint of the original article published by the 
Institute of Mathematical Statistics in The Annals of Probability, 
2011, Vol. 39, No. 5, 1983-2017. This reprint differs from the original in 
pagination and typographic detail. 



1 



2 



R. KENYON 



(and therefore, each component is a tree plus an edge: a unicycle). Here the 
weight of a configuration is the product of 2 minus the trace of the holonomy 
around each of the cycles. 

By varying the connection on the vector bundle, one constructs in this 
way many natural measures on CRSFs. If the monodromy of the underlying 
connection is unitary, then these measures are determinantal processes for 
the edges (or g-determinantal in the case of an SL2C-connection). 

We give here a number of applications of these results. Further applica- 
tions can be found in the papers [6, 9] and [13]. 

1.1. Spanning trees. One of the main applications is to the study of 
spanning trees, in particular spanning trees on planar graphs or graphs on 
surfaces. By a judicious choice of vector bundle connection, the natural 
probability measure on CRSFs can be made to model the uniform spanning 
tree measure conditioned on having certain "boundary connections." As 
simple examples we compute the probability that, on a finite planar graph, 
the branch of the uniform spanning tree connecting two boundary points zi 
and Z2 passes left of a given face or a given two faces. These probabilities 
are given in terms of the Green's function for the standard Laplacian on the 
graph. 

1.2. -functions. Suppose ^ is a graph embedded on an annulus. On such 
a graph, a 6 -function is a multi- valued harmonic function whose "analytic 
continuation" around the annulus is a constant times the original function. 
In other terminology (defined below), it is a harmonic section of a flat line 
bundle. The constant is called the multiplier of the 9 function. We prove 
that if k is the largest integer such that one can simultaneously embed k 
pairwise disjoint cycles in Q, each winding once around the annulus, then Q 
has at most 2k — 1 ^-functions. These ^-functions have distinct positive real 
multipliers; these multipliers are closed under inverses (if A is a multiplier 
then so is A~^). For generic conductances, G will have exactly 2k — 1 9- 
functions; we conjecture that this is always the case. The multipliers have 
a probabilistic meaning for CRSFs: see below. 

1.3. Graphs on surfaces. For a graph embedded on a surface (in such 
a way that complementary components are contractible or peripheral an- 
nuli), a natural probability model is the uniform random CRSF whose cycles 
are topologically nontrivial (not null-homotopic) . Such CRSFs are called in- 
compressible. The cycles in a CRSF give a finite lamination of the surface, 
that is, a finite collection of disjoint simple closed curves. The Laplacian 
determinant on a flat line bundle on the surface counts CRSFs weighted by 
a function of the monodromy of the connection. By considering the Lapla- 
cian determinant as a function on the representation variety (consisting of 



SPANNING FORESTS AND THE VECTOR BUNDLE LAPLACIAN 



3 



flat connections modulo conjugacy), one can extract the terms for each pos- 
sible topological type of finite lamination. In particular, one can study the 
uniform measure on incompressible CRSFs. For example, the number of 
components for a uniform incompressible CRSF on a graph on an annulus 
is distributed as a sum of a finite number of independent Bernoulli ran- 
dom variables, with biases given by where A runs over the ^-function 
multipliers defined above. See Corollary 1 below. 

Given a Riemann surface one can take a sequence of finer and finer graphs 
adapted to the metric so that the potential theory on the graph and on the 
Riemann surface agree in the limit. In this case, one can show [9] that the 
uniform incompressible CRSF has a scaling limit, whose distribution only 
depends on the conformal structure of the underlying surface. This is similar 
to the theorem of Lawler, Schramm and Werner on conformal invariance of 
the uniform spanning tree [16]. We compute here for the annulus and the 
square torus the exact distribution of the number and homology class of the 
cycles of an incompressible CRSF in the scaling limit. 

1.4. Monotone lattice paths. The papers [7] and [14] studied CRSFs on 
a north/east-directed m x n grid on a torus. The distribution of the num- 
ber of cycles was shown to have a highly nontrivial structure as a function 
of m,n. Here the line bundle Laplacian gives quantitative information about 
the model which was not available in [7]. In particular, we obtain an exact 
expression for the generating function for the number and homology type 
of the cycles. For the nx n torus, for example, we show that the number of 
cycles tends to a Gaussian as n — )■ oo with expectation and variance 



2. Background. The uniform probability measure on spanning trees on 
a graph, called the UST measure, has for the past 20 years been a remarkably 
successful and rich area of study in probability theory. Pemantle [19] showed 
that the unique path between two points in a uniform spanning tree has the 
same distribution as the loop-erased random walk between those two points.^ 
Wilson [20] extended this to give a simple method of sampling a uniform 
spanning tree in any graph. Burton and Pemantle [3] proved that the edges 
of the uniform spanning tree form a determinantal process (see definition 
below). This allows computation of multi-edge probabilities in terms of the 
Green's function. 

The UST on 7? received particular attention due to the conformal in- 
variance properties of its scaling limit. Pemantle showed that almost surely 



^The loop-erased random walk (LERW) between a and b is defined as follows: draw 
the trace a simple random walk from a stopped when it reaches b, and then erase from 
the trace all loops in chronological order. What remains is a simple path from a to b. 
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the UST in 1? has one component. In [10], we proved that the expected 
length of the LERW in Z^, and therefore a branch of the UST of diameter n 
was of order n^^^, a result predicted earher by conformal field theory [17]. 
In [12], we showed that the scaling limit of the "winding field" (describing 
how the branches of the UST wind around faces) was a Gaussian free field. 
Further conformally invariant properties were proved in [11]. In [16], Lawler, 
Schramm and Werner proved that the Z^-LERW converges to SLE2, and the 
peano curve winding around the Z^-UST converges to SLEg. 

As a result of these works, we have a decent understanding of the scaling 
limit of the UST on 1? . However, some important questions remain. For 
example, how are different points in the UST connected: given a set of 
points in the plane, what is the topology of their tree convex hull, that 
is, the union of the branches of the UST connecting them in pairs? Can 
one compute various connection probabilities, for example, the probability 
that the LERW from (0, 0) passes through the points t^i, ^2, . . . , t^n in order? 
What is the distribution of the "bush size," the finite component of the 
tree obtained by removing a single random edge? While many of these can 
be answered in principle using SLE techniques, in practice one must solve 
a hard PDE. 

Many of these questions can be answered with the bundle Laplacian. Some 
of these are illustrated below. We will discuss how these results can be used 
to compute various connection probabilities for the UST in [13], and discuss 
the connection with integrable systems in [6]. 

3. Vector bundles on graphs. 

3.1. Definitions. Let Q he a finite graph. Given a fixed vector space V, 
a V -bundle, or simply a vector bundle on Q is the choice of a vector space Vy 
isomorphic to V for every vertex v of Q. A vector bundle can be identified 
with the vector space Vg := ^^V^ = yl^L A section of a vector bundle is 
an element of Vg. 

A connection $ on a ^-bundle is the choice for each edge e = vv' of G of an 
isomorphism (pyy/ between the corresponding vector spaces (j)^^' : Vy — )• V^i , 
with the property that (f)^^' = (f>~,\;. This isomorphism is called the parallel 
transport of vectors in Vy to vectors in Vy/. Two connections <I>, are said to 
be gauge equivalent if there is for each vertex an isomorphism ijjy '.Vy ^ Vy 
such that the diagram 

it , 

Vv — ^ K' 

commutes. In other words ^' is just a base change of <I>. Given an oriented 
cycle 7 in Q starting at the monodromy of the connection is the ele- 
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ment of End(Vu) which is the product of the parallel transports around 7. 
Monodromies starting at different vertices on 7 are conjugate, as are mon- 
odromies of gauge-equivalent connections. 

A line bundle is a F-bundle where V = C, the one-dimensional complex 
vector space. In this case, if we choose a basis for each C then the parallel 
transport is just multiplication by an element of C* = C \ {0}. The mon- 
odromy of a cycle is in C* and does not depend on the starting vertex of 
the cycle (or gauge). 

3.2. The Laplacian. The Laplacian A on a ^-bundle with connection <I> 
is the linear operator A:Vg ^ Vg defined by 

Af{v) = Y,fiv)-<Pv'vf{v'), 

where the sum is over neighbors v' of v. 

If we assign to each edge a positive real weight (a conductance) Cy^' = c^'y, 
the weighted Laplacian is defined by 

Af{v)=Y,Cyy'{f{v)-(l)y>yf{v')). 
v' 

Note that if the vector bundle is trivial, in the sense that (j)yyi is the 
identity for all edges, this is the classical notion of graph Laplacian (or more 
precisely, the direct sum of dimT/ copies of the Laplacian). 

Here is an example. Let Q = K3 with vertices {vi,V2,V3} . Let ^> be the 
line bundle connection with (py-y- = Zij € C*. Then in the natural basis, A 
has matrix 

/ 2 -Zi2 -zisX 
(1) A=[-Z^^ 2_^ -Z23 . 

3.3. Edge bundle. One can extend the definition of a vector bundle to 
the edges of ^. In this case, there is a vector space Ve = V for each edge e 
as well as each vertex. One defines connection isomorphisms cpye = (p'^y for 
a vertex v and edge e containing that vertex, in such a way that if e = vv' 
then <j)yyf = (pey' ° 4>ye, where (pyy' is the connection on the vertex bundle. 

The vertex/edge bundle can be identified with yl^^l+l-^l = Vg® Ve, whe- 
re Ve is the direct sum of the edge vector spaces. 

A 1-form (or cochain) is a function on oriented edges which is antisym- 
metric under changing orientation. If we fix an orientation for each edge, 
a 1-form is a section of the edge bundle, that is, an element of yl^L We 
denote by A^(^,<I>) the space of 1-forms and A^{Q,^) the space of 0-forms, 
that is, sections of the vertex bundle. 
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We define a map d:A^{g,<^) ^ A\g,^) by df{e) = (j^yefiy) - ^^ef{x) 
where e = xy is an oriented edge from vertex x to vertex y. We also define 
an operator d* : — )■ A'' as follows: 

d*u{v) = ^ 0ei.w(e), 

e=v'v 

where the sum is over edges containing v and oriented toward v. Despite the 
notation, this operator d* is not a standard adjoint of d unless (j)ev and (p^e 
are adjoints themselves, that is, if parallel transports are unitary operators 
(see below). 

The Laplacian A on A*' can then be defined as the operator A = d*d: 
d*df{v)= ^evdfie) 

e=v'v 

= X] 'Pev{4>vef{v) - (l)v'ef{v')) 
e=v'v 

v' 

= Afiv). 

We can see from the example (1) above on that A is not necessarily 
self-adjoint. However, if (f)„,i is unitary: = cp*,^ then d* will be the adjoint 
of d for the standard Hermitian inner products on V^^^ and V^^^, and so in 
this case A is a Hermitian, positive semidefinite operator. In particular on 
a line bundle if | | = 1 for all edges e = vv' , then A is Hermitian and 
positive semidefinite. 

4. Spanning forests associated to a line bundle. 

4.1. Cycle-rooted spanning forests (CRSFs). For a line bundle {Q,^) on 
a finite graph we have a combinatorial interpretation of the determinant 
of the Laplacian in terms of cycle-rooted spanning forests. A cycle-rooted 
spanning forest (CRSF) is a subset S of the edges of Q, spanning all vertices 
(in the sense that every vertex is the endpoint of some edge) and with the 
property that each connected component of S has as many vertices as edges 
(and so has a unique cycle). See Figure 1. 

An oriented CRSF is a CRSF in which we orient each edge in such a way 
that cycles are oriented coherently and we orient the edges not in a cycle 
toward the cycle. An oriented CRSF is the same as a nonzero vector field 
on which is the choice of a single outgoing edge from each vertex. 

A cycle-rooted tree (CRT), also called unicycle, is a component of a CRSF 
(oriented or not). 
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Fig. 1. A CRSF on a square grid. 



4.2. The matrix-CRSF theorem. The following theorem is due to For- 
man [5]. 

Theorem 1 [5]. For a line bundle on a connected finite graph, 
detA= '[[{2-w-l/w), 

CRSFs cycles 

where the sum is over all unoriented CRSFs C , the product is over the cycles 
of C , andw,l/w are the monodromies of the two orientations of the cycle. 

Recall that for a line bundle the monodromy of an oriented cycle does 
not depend on the starting vertex. Note that we could as well have written 
det A = X^ocRSFs ncyciesC-*^ ~ where now the sum is over oriented CRSFs 
(OCRSFs). 

Forman uses an explicit expansion of the determinant as a sum over the 
symmetric group, and a careful rearrangement of the terms. We give a differ- 
ent proof using Cauchy-Binet formula which allows us to bypass this step. 

Proof of Theorem 1 . We use the Cauchy-Binet formula, det A = 
detd*d = det{B* B) , where B is a maximal minor of d (i.e., if d is n x m 
then B runs over all n x n submatrices of d obtained by choosing n columns). 

Columns of B index vertices and rows of B index edges of Q. Since B is 
square, it corresponds to a selection of n edges (where n is the number of 
vertices of ^). 

If the set of edges in B contains a component C with no cycles, say of size 
\C\ = k, this component necessarily has k — 1 edges. So B maps a /c-dimen- 
sional subspace of the vertex space to a A; — 1-dimensional subspace of the 
edge space and so deti? = 0. Therefore to have a nonzero contribution to 
the sum every component of B has at least one cycle, and, since the total 
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number of vertices equals the total number of edges, every component must 
have exactly one cycle. We have proved that the sum is over CRSFs. 

If B has only one component which is a CRT, orient its edges so that all 
edges are oriented toward the unique cycle, and along the cycle choose one 
of the two coherent orientations. Write 

detB= ^ sgn(cr)Si„(i) 

The only nonzero terms in this expansion are those in which vertex i is 
adjacent to edge cr{i), that is, the two terms which are consistent with one 
of these two orientations. Thus 




where the products over j are over edges not in the cycle (called bush edges) , 
and those over i are for the two cyclic orientations (indices are taken cycli- 
cally). If the cycle has odd length, sgn((Ti) = sgn((T2) and the sign change 
in front of the second product in (2) is due to switching the sign on the 
odd number of edges of the cycle. If the cycle has even length, the sign 
change is due to sgn(f7i) = — sgn((T2) (a cyclic permutation of even length 
has signature —1). 
We similarly have 

deti?* = sgn(ai)(-l)" J] ^-.-.(U " 11 '^^^-m)' 

jSbushes iScycle iScycle 

and multiplying these two and using (pev4^ve — 1; 

.et..B.(.,-.,(i.-^)..-ii-^. 

where ^ = rijecycie 'PviVi+i which is the monodromy of the cycle. 

If each component of is a CRT, then det -B*-B is the product over 
components of the above contributions. □ 

5. Measures on CRSFs. In the case of a line bundle $ with unitary 
connection. Theorem 1 provides a definition of a natural probability mea- 
sure /i4> on CRSFs. A CRSF has probability equal to ^11(2 ~ ^ ~ where 
the product is over its cycles. The constant of proportionality Z, or parti- 
tion function, is then just det A. Note that as long as the (unitary) bundle 
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is not trivializable (gauge equivalent to the trivial bundle) then some cycle 
has nontrivial monodromy and so det A > 0. 

The goal of this section is to show that the measure is determinantal for 
the edges. 

5.1. Harmonic sections. The Laplacian is the composition: 

If A has full rank, then A^(t/, ^) splits as 

=Im((i)eKer(d*). 

A section f £ is harmonic if A/ = 0. If A does not have full rank, 
there are nontrivial harmonic sections and the above is no longer a direct 
sum. A 1-form to G is harmonic if it is in the intersection lm{d) n Ker(d*). 
Such a form is both exact (lo = df for some /) and co-closed [d*uj = 0). 

Proposition 1. When A has full rank, the projection along Ker(d*) 
from A^ to the space of exact \ - forms \m.{d) is given by the operator P = dgd* 
where g = A~^ is the Green's function for A on A''. 

Proof. Note that P = dgd* is zero on co-closed forms and the identity 
on exact forms: P{df) = dgd*df = dgAf = df. □ 

This projection operator P plays a role below. 

5.2. Determinantal processes. Let ;U be a probability measure on = 
{0, 1}"". For a set of indices C {1, 2, . . . ,n}, define the set 

Es = {{xi, . . . ,x„) G r2 1 = 1 Vs G 5}. 

We say ^ is a determinantal probability measure (see, e.g., [8]) if there is an 
nxn matrix K with the property that for any set of indices S C {1, 2, . . . , n}, 
we have 

fi{Es) = det[{Kij)ij(zs]. 

That is, principal minors of K determine the probability of events of type Es- 
The matrix K is called the kernel of the measure. A simple example of 
a determinantal process is the product measure, in which case K is a di- 
agonal matrix with diagonal entries in [0,1]. Another well-known example 
of a determinantal process is the uniform spanning tree measure on a finite 
graph [3]. Here {1,2,..., n} index the edges of a connected graph. The ker- 
nel K is the so-called transfer- current matrix, defined hy K = dGd* , where G 
is the Green's function, see [3]. 

A simple inclusion-exclusion argument shows that individual point prob- 
abilities are also given by determinants: 
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Proposition 2. Let X = {xi,X2,-.-,Xn) Then 
M^) = (-ir-l^ldet(diag(X)-if), 

where diag(X) is the diagonal matrix whose diagonal entries are 1 — Xi, and 
1^1 =E^=l^^• 
5.3. A determinantal process on CRSFs. We will need the following well- 
known lemma. 

Lemma 1. If ^) is a block matrix with A, D square and D invertible 
then 



A B 
C D 



det 

Proof. This follows from 



^ = det D det{A - BD-^C). 



A B\_(A-BD^^C BD-^\ f I 

c Dj ~y I J yc D J ■ □ 

Recall the definition of the projection P from 1-forms to exact 1-forms 
(Proposition 1). 

Theorem 2. For a line bundle $ with unitary connection, with respect 
to the measure /i$ the edges of the CRSF form a determinantal process with 
kernel P = dgd* . 

Proof. Let ^ be a graph with n vertices and m edges. Let {ei , . . . , e„,} C 
E be the edges of a CRSF 7. Write the matrix for d and d* so that the first n 
edges are ei, . . . , e^. Then d = (^^) where di is n x n and d2 is (m — n) x n. 

From Proposition 2, the probability of 7 is Pr(7) = (—1)"*"" det{X — dgd*) 
which by Lemma 1 with A = X,B = d,D = A and C = d* can be written 



Pr(7) = (-ir 



di 
' '* ' ' detd^di 



det A det A 

By Theorem 1 this is exactly the probability of 7. □ 

Among other things, this theorem along with Proposition 2 allows us to do 
exact sampling from the measure as follows (see [8]). Pick an edge ei; it 
is present with probability P(ei, ei). Take another edge 62; if ei is present, 62 
will be present with probability 

Pr(e2|eO = T^^detf^J^^'^i) ^f^^'^^S 
P(ei,ei) \P{e2,ei) P(e2,e2j 
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If ei is not present, 62 will be present with probability 




P(ei,ei) P(ei,e2) 
P{e2,ei) P{e2,e2) 



) 



) 



and so on. 

5.4. Small nionodromies. An important limit, or actually set of limits 
of these determinantal processes, is when the monodromies tend to 1. Even 
though det A tends to zero, the probability measures may converge. We 
get determinantal processes supported on CRSFs on the graph with trivial 
line bundle. 

Theorem 3. Choose an orientation for each edge and suppose the paral- 
lel transport is (j)^ . = 6***^^ for the jth edge in the direction of its orientation. 
Fix the Cj and let t ^ 0. In the limit t — )• the determinantal process yU^j 
tends to a determinantal process hq = ^o({cj}) supported on CRSFs with 
a single component, that is, CRTs. The probability of a CRT is proportional 
to (X^Cj)^, where the sum is over oriented edges in the unique cycle. 

Proof. The monodromy of a loop is -u; = e**^'^^ . We have 2 — w — ^ = 
2 — 2cos(t ^Cj) = Cj)^ + 0{t^). In the limit t — )■ 0, the partition function 
is Z = 0{t^) so only CRSFs with one component contribute. □ 

6. An application io Q X "Ln- Let ^ be a finite graph with m vertices 
and Z„ the n-cycle. Let be the product graph whose vertices are ^ x Z„ 
and edges connect (x, j) to + 1 mod n) and {x,j) to {x' ,j) when x,x' 
are neighbors in Q. 

We compare the uniform spanning tree on Hn and the uniform cycle- 
rooted tree whose cycle winds around Z„. The minimum cut set of a span- 
ning tree is the set of edges not in the tree such that when added to the tree 
make a cycle winding nontrivially around Z„. 

Fix z G C with \z\ = 1. On the product graph Hn let = 1 except 
when V = {x,n— 1) and v' = {x, 0) in which case (j)yyi = z. Then cycles have 
trivial monodromy unless they wind nontrivially around the Z„ direction. 
So a CRSF has nonzero probability only if all of its cycles wind nontrivially 
around Z„ (possibly many times). 

Because of the product structure the eigenvectors of the Laplacian on Hn 
are products of eigenvectors of the standard Laplacian on Q and the line 
bundle Laplacian on ^7^. Thus, the eigenvalues of the Laplacian on Hn are 
sums of eigenvalues of Ag (the standard Laplacian) and Az„ (the line bundle 
Laplacian). The eigenvalues of Az„ are 2 — — (^~^ where C ranges over the 
roots of = z (the corresponding eigenvectors are exponential functions). 
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The Laplacian determinant of Hn is then 

detA= n det(Ae + 2-C-r^)= [] H^^ + ^-C-C"') 

(3) 

= (2-^-^-1) J] JJ(A + 2-C-rt, 

where A runs over the eigenvalues of the Laplacian on Q. 

Compare this to det' Aq, the product of nonzero eigenvalues of the stan- 
dard Laplacian Aq on Hn- 

det'Ao= W A + 2-C-r^ 

(C,A)^(1,0) 

(4) = n 2-c-c-i n 

^'1=1 XjLQ 



Theorem 4. As n ^ oo a fi^-random CRSF has one component with 
probability tending to 1, that is, is a CRT. Its uniQue cycle winds avound 
once with probability tending to 1. The ratio Rn of the number of such 
CRTs and number of spcLTiTiifiQ trees satisfies lini^^Qo — — . The expected 
length IE(L„) of the cycle in such a random CRT of Hn satisfies 

lim ^^-^"^ = 1 

n^oo nE{Sn) 

where is the expected size of the minimum cut set of a random span- 

ning tree of Hn . 

Proof. The first two statements can be seen as follows. Take n large 
and consider the part of the configuration in a piece Q x [i,i + N] for large A^. 
We claim that the number of spanning tree configurations restricted to this 
subgraph, as a function of A^, has a strictly larger exponential growth rate 
than that of the number of spanning forest configurations in this subgraph 
with two or more components ( "strands" ) each of which intersects both ends 
X {i} and G x {i + N}. The growth rate of spanning trees can be com- 
puted as the leading eigenvalue of nonnegative finite matrix T, the transfer 
matrix.^ It is easily seen that the transfer matrix T is primitive, that is, has 



''The states of the transfer matrix are all possible forests (edge subsets) oi Q x {i} 
occurring in a spanning tree oi Q x Z along with a partition of the components in this 
forest, which describes which components of this forest are connected to each other to the 
left, that is, in the part of the tree in Q x ( — oo,i — 1]. 
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the property that some power is strictly positive. On the other hand, the 
transfer matrix for spanning forest configurations with two or more strands 
is a matrix which can be obtained from T by setting some of its positive 
entries to (those entries whose components are ah connected to the left). 
This strictly decreases its leading eigenvalue. This proves the claim and the 
first two statements of the theorem. 

If we divide (3) by 2 — z — and let z tend to 1 the expression counts 
CRSFs with one component, that is, CRTs, with a weight A;^ if they wind k 
times around Z/nZ, because 

2 2 ^ 

lim r- = . 

z^i 2 - z- z-^ 

As n gets large CRTs which wind more than once around have probability 
tending to zero. Line (4) is the number of rooted trees; dividing by nm, the 
number of locations for the root, gives the number of trees, which by (3) and 
the above remarks is ^ times the number of CRTs winding once around, 
plus errors tending to zero as n — )■ oo. This proves the third statement. 

The following two sets A and B are in bijection: A is the set of CRTs 
whose cycle winds once around along with the choice of an edge on this 
cycle. B is the set of spanning trees with a choice of a complementary edge 
in the minimum cut set. The bijection consists in adding the edge to the 
tree. Therefore, the expected length of the unique cycle in a random CRT is 

\A\ _ \B\ _ \B\ |trees| 

|CRTs| ~ |CRTs| ~ |trees| |CRTs|' 

By the above, 1^^^^. □ 

7. Weighted and/or directed graphs. 

7.1. Edge weights. Putting weights on the edges is a minor general- 
ization; let c:E ^ ]R>o be a conductance associated to each edge, with 
Cvv' = Cy'v We then define 

^f{v) = '^Cyyi{f{v) - (l)v'vf{w)). 

In other words, letting C be the diagonal matrix, indexed by the edges, 
whose diagonal entries are the conductances, we have 

Theorems. ^'sk,[^ = Y^(j^^-^^\{^C{e)\{(2-w-^) where the first prod- 
uct is over all edges of the configuration and the second is over cycles, and w 
is the monodromy of the cycle as before. 

This follows directly from the proof of Theorem 1 above. 
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7.2. Directed graphs. We can also have weights which depend on direc- 
tion, so that each edge has two weights, one associated to each direction. 
This makes G into a Markov chain, in which the transition probabihties are 
proportional to the edge weights. 

In this case the Laplacian has the same form 

^f{v) = ^ - (pv'vfiv')), 

v' 

except that Cw' 7^ c^'v in general. The operator d : A*^ — )■ can be defined 
by 

dfie) = (pvefiv) - (pv'efiv') 

as in the unweighted case, where e = v'v, and d* by 

d*{uj){v)= ^ c„„'(^et,w(e), 

e=v'v 

where the sum is over edges e = v'v containing v. 

In this case, each oriented CRSF has a different weight. Forman's theorem 
in this case is the following. 

Theorem 6 [5]. 

detA= Yl n n ^^(7)(l-^(7)), 

OCRSFs eSbushcs cycles 7 

where the sum is over oriented CRSFs, the first product is over the edges in 
the hushes (i.e., not in the cycles), oriented toward the cycle, and the second 
product is over oriented cycles, C{'j) is the product of semiconductances 
along 7 and w is the monodromy of 7 . 

An application is in Section 13 below. 

8. Laplacian with Dirichlet boundary. Let ^ be a graph with line bundle 
and connection $ and S a subset of its vertices, which play the role of 
boundary vertices. We can define a Laplacian Ag^5 with Dirichlet boundary 
conditions at S as follows: for / S V^\^ and f S ^ \ S, 

^g,sf{v) = {degv)f{v)- Yl <t>v'J{v'). 

v' S:Q\S ,v' 

This is of course just the Laplacian Ag restricted to the subspace V^"^^ 
of functions which are zero on S and projected back to this subspace. As 
a matrix, it is just a submatrix of Ag. 

An essential CRSF of (G, S) is a set of edges in which each component is 
either a tree containing a unique vertex in 5 or a cycle-rooted tree containing 
no vertices in 5. 
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Theorem 7. det Ag.^ is the weighted sum of essential CRSFs, where 
each configuration has weight ncRTs(^ — w — w~^), where w is the mon- 
odromy of the cycle. 

This generalizes the matrix-tree theorem^ (when $ is the identity and 
l^l = 1) and the matrix-CRSF theorem (when IS"! = 0). 

Proof of Theorem 7. The proof follows the same lines as the proof 
of Theorem 1. Now d is an m x n matrix with n = \Q\ — \S\. A maximal 
minor B of d corresponds to a choice of n edges. Each component of B 
is either a CRT, in which case its weight is computed as before, or a tree 
connecting some vertices of Q\S with a single vertex of S. In this case, the 
determinant of B (we mean, that part of B coming from this component) 
has a single nonzero term in its expansion, and the corresponding term in B* 
is its inverse. So each tree component counts 1. □ 

Again when = 1 there is a determinantal measure associated with A. 

9. Bundles with SL2C connections. We can extend many of the above 
results to the case of a C^-bundle with SL2C-connection. 

9.1. Q- determinants. Given A £ GL2(C) define A = (det^)^""^, the ad- 
jugate of A. That is, if A = (" ^) then i = (^ ^). Note that 

(5) A + A = Tr{A)I 
is a scalar. 

Let M be a matrix with entries in GL2(C). M is said to be self-dual if 
Mij = Mji. In particular, diagonal entries must be scalar matrices. 

For a matrix M with noncommuting entries, there are many possible ways 
one might define its determinant: however these all involve the expansion 

(6) detM = ^ sgn(a)[Mi,(i) • • • M„,(„)]^, 

where [Mi„(^i^ ■ ■ ■ Mna(n)]cT denotes the product of the terms Mjo.(j) after they 
have been rearranged in some specified order depending on a. Each of these 
"order functions" determines a different possible determinant. 

In the current case, there is one very natural condition to put on these 
orders. If a is written as a product of disjoint cycles, then in the corre- 
sponding order the terms in each cycle of a should appear consecutively. 



^One version of the classical matrix-tree theorem is that the determinant of the matrix 
obtained by removing a single row and column v from A is equal to the number of spanning 
trees rooted at v. 
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Moreover, if o", a' are two permutations with the same cycles, except that 
some of the cycles are reversed, then the corresponding order should have 
the cycles appear in the same relative order. 

The advantage of this for self-dual matrices is that the product of entries 
along a cycle is the adjugate of the product along the reversed cycle; so that 
the sum of these is a scalar by (5). Thus by grouping permutations into 
sets with the same cycles up to reversals, one arrives at a product of scalar 
matrices. The sum (6) with appropriately rearranged products then yields 
a 2 X 2 scalar matrix q ■ I, and the number q is defined to be the determinant. 

The Q-determinant of self-dual matrix M [18] is defined in precisely this 
way. We define 

Qdet(M)= sgn(^) J] ^^^H' 

cr^Sn cycles 

where the sum is over the symmetric group, each permutation a is written 
as a product of disjoint cycles, and tic{w) is the trace of the product of the 
matrix entries in that cycle. If we group together terms above with the same 
cycles — up to the order of traversal of each cycle — then the contribution 
from each of these terms is identical: reversing the orientation of a cycle 
does not change its trace. So we can write 

c 

Qdet{M) = Y,{-^r^''ll^r{wi), 

1=1 

where the sum is over cycle decompositions of the indices, c is the number of 
cycles, and Wi is the monodromy (in one direction or the other) of each cycle. 
Here tr is equal to the trace for cycles of length at least 3; cycles of length 1 
or 2 are their own reversals so we define tr(u;) = ^tT{w) for these cycles. 
As an example, let A = al,C = cl and i3 = (^J ). Then 

Qdet(^^ c) = tr(^)to(C)-G(SS) = ac- (6164- 62fe3)- 

Note that if M is a self-dual nx n matrix then ZM, considered as a 2n x 
2n matrix is antisymmetric, where Z is the matrix with diagonal blocks 
(^-^ 0) 'iiid zeros elsewhere. The following theorem allows us to compute 
Q-determinants explicitly. 

Theorem 8 [18]. Let M be an n x n self-dual matrix with entries in 
GL2(C) and M' the associated 2n x 2n matrix, obtained by replacing each 
entry with the 2 x 2 block of its entries. Then Qdet(M) = Pf(ZM'), the 
Pfaffian of the antisymmetric matrix ZM' . 

9.2. Laplacian determinant. For a C^-bundle on a graph Q with SL2C 
connection A is a self-dual operator (in the above sense). 
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Theorem 9. For a C'^ -bundle on a connected graph Q with SL2C con- 
nection, 

QdetA= ^ JJ 2-tr(u;), 

CRSFs cycles 

where tx{w) is the trace of the monodromy of the cycle. 

Proof. The proof generalizes the proof of Theorem 1. We use as before 
the Cauchy-Binet formula which holds for general Q-determinants (Theo- 
rem 10 below). 

Suppose G has n vertices. Write Qdet A = Qdet{d*d) = Y^^QdetB*B. 
Here B runs over choices of n edges of ^. If -B corresponds to a set of edges 
which has a component with no cycle, that is, a component with k — 1 edges 
and k vertices, the corresponding matrix B is singular (in fact has a block 
form with nonsquare blocks) and so Qdet B = V det B = 0. 

Thus, as before, each component of B must have exactly one cycle. We 
now work directly with B*B. Note that this is the Laplacian on the subgraph 
defined by B. Suppose that B has only one component; the general case is 
similar. A nonzero term in the expansion of the determinant corresponds to 
a decomposition of a into cycles; the only such terms which are nonzero are 
when adjacent vertices are paired, each vertex is fixed, or the vertices on the 
unique cycle are advanced or retreated in the direction of that cycle and the 
other vertices are fixed. 

Consider the case when some vertices are paired; if i,j are paired in B*B 
then B maps i and j to the same edge Cij, so this term does not contribute. If 
no terms are paired, the expansion of the Q-determinant for B*B therefore 
has exactly four terms. It has two terms in which a fixes the bush vertices and 
advances the cycle vertices around the cycle in one of the two orientations. 
It also has two terms in which all vertices are fixed: the vertices in the cycle 
advance to the edges in one of the two directions and then retreat back to 
their original positions. The signs work out as before. □ 

Theorem 10 (Noncommutative Cauchy-Binet formula). 

Qdet(MiM2) = J]] Qdet(BiB2), 

where Bi runs over maximal minors of Mi, and B2 is the corresponding 
minor of M2 ■ 

Proof (Sketch). Suppose Mi is n x m with m> n. The standard proof 
only relies on the multilinearity of the determinant. Write columns of M1M2 
as linear combinations of columns of Mi with coefficients in M2 (acting on 
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the right). Use multihnearity to write this as a sum of m" Q-determinants; 
the Q-determinant of one of these is zero if two of the same columns are 
used. Group the remaining m choose n nonzero terms into products 
QdetSiQdet^s. □ 

10. Graphs on surfaces. 

10.1. Definitions. Let S be an oriented surface, possibly with boundary, 
and Q a graph embedded on T, in such a way that complementary compo- 
nents are contractible or peripheral annuli (i.e., an annular neighborhood of 
a boundary component). We call the pair {Q,T,) a surface graph. 

For each boundary component there is a "peripheral" cycle on Q, consist- 
ing of those vertices and edges bounding the same complementary compo- 
nent as the boundary component. 

An finite lamination or simple closed curve system on S is (the isotopy 
class of) a finite set of pairwise disjoint simple closed curves on S each of 
which has nontrivial homotopy type (i.e., no curve bounds a disk). We allow 
two or more curves in the lamination to be isotopic. 

10.2. Bundles on surface graphs. Given a surface graph {G,T,) a vector 
bundle on Q with connection $ is flat if it has trivial monodromy around 
any face (contractible complementary component of S; peripheral cycles do 
not bound faces). In this case, given a closed loop starting at a base point v, 
the monodromy around a loop only depends on the homotopy class of the 
(pointed) loop in 7ri(S), and so the monodromy determines a representation 
of 7ri(S) into End{V). This representation depends on the base point for tti; 
choosing a different base point will conjugate the representation. 

Conversely, let p G Hom(7ri(S), End(y)) be a representation of '/ri(S) into 
End(y); for a fixed base point there is a unique flat line bundle, up to 
gauge equivalence, with monodomy p. It is easy to construct: start with 
trivial parallel transport on the edges of a spanning tree of Q. Now on 
each additional edge, adding it to the tree makes a unique cycle; define 
the parallel transport along this edge to be the image of the homotopy class 
of this cycle under p. Thus, we have a correspondence between flat bundles 
modulo gauge equivalence and a space X = IIom(7ri(5]),End(y))/End(y), 
the space of homomorphisms of vri(S) into End(y) modulo conjugation. This 
space X has the structure of an algebraic variety: it can be represented using 
variables for the matrix entries of a set of generators of 7ri(S), modulo an 
ideal corresponding to the relations in ■7ri(5]) and the conjugation relations. 
The variety X is called the representation variety of vri(S) in End(y). 

When dealing with flat bundles it is natural to restrict our attention to 
CRSFs each of whose cycles is topologically nontrivial, since other CRSFs 
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have zero weight in det A. We call these incompressible CRSFs. The cycles 
in an incompressible CRSF form a finite lamination. 

We now restrict to ^ = and SU2C-connections. For a flat bundle, det A 
is a regular function on the representation variety X = Hom(7ri(S), SU2C)/ 
SU2C (i.e., it is a polynomial function of the matrix entries). Remarkably, by 
varying the representation it is possible to extract from det A those terms for 
incompressible CRSFs whose cycles have any given set of homotopy types. 
The following theorem is due to Fock and Goncharov [4]: 

Theorem 11 ([4], Theorem 12.3). // S has nonempty boundary, the 
products 

n Tr(7) 

cycles 7 

over all finite laminations form a basis for the vector space of regular func- 
tions on the representation variety Hom(7ri(S), SU2C)/SU2C. 

In our case, det A is a regular function on X; hence it can be written 
detA = J^^CL W TrH, 

L cycles 

where G Z, the sum is over finite laminations L, the product is over 
curves 7 in L and w is the monodromy of 7. By an integral change of basis 
(in fact an upper triangular integer matrix, if we order isotopy types of finite 
laminations by inclusion) we can write this as 

detA = ^4 \{{2-Tt{w)), 

L cycles 

and therefore by Theorem 9, is the desired weighted sum of CRSFs with 
lamination type L. 

The actual extraction of a coefficient can be done as follows. There is a na- 
tural measure /i on X, essentially just the product of k copies of Haar 
measure on SU2C, where k is the rank of the free group 7ri(S). The reg- 
ular functions discussed above form a basis for the Hilbert space L'^{X,fi) 
(not an orthonormal basis but obtained from an orthonormal basis by an 
upper-triangular linear transformation). Thus, one can obtain cl above by 
an integration of det A against a particular function over X. 

10.3. Examples. 

10.3.1. Annulus. The simplest example of a surface graph, after a graph 
on a disk, is a graph on an annulus. 
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Suppose is a surface graph on an annulus. Since vri(S) = Z is 

Abelian, it will suffice to work with a line bundle (a Z subgroup of SL2C 
will have the same traces as a Z subgroup of diagonal matrices in SL2C, 
which is equivalent to taking a line bundle connection). Let z €C* be the 
monodromy of a flat line bundle. 

Then P{z) = det A is a Laurent polynomial in z. It is reciprocal: P{z) = 
P{l/z) by Theorem 1. 

Theorem 12. P{z) is a reciprocal polynomial with real and positive 
roots and a double root at z = 1. 

We emphasize here that the roots of P are unrelated to the eigenvalues 
of A. Rather they are special values of the monodromy z for which A is 
singular. We conjecture that the roots are distinct. 

Proof of Theorem 12. When z = l, the line bundle is trivializable 
and A is the usual graph Laplacian. Thus, P{1) = 0. However, there exists 
a CRSF on Q with a single component winding once around S: just take any 
cycle of this type and complete it to a CRT. So in Theorem 1, the coefficient 
of {2 — z — l/zY is nonzero, and z = 1 is a double root (i.e., not of higher 
order) . 

We prove reality by a deformation argument. At each real root r > 1 of P, 
A has a kernel Wr- We claim that dim = 1. If it were of larger dimension, 
Ist /i, /2 S Wr be independent. These are real and harmonic for A, and so lift 
to actual harmonic functions on the embedded graph (^, S) on the universal 
cover S of S. These functions have the property that f{x-\-l) = rf{x) where 
"+1" represents the deck transformation, that is, on a path winding once 
around the annulus the values of /i, /2 are multiplied by r. 

Let V be a vertex on the boundary of T, (i.e., on a peripheral cycle) and /s 
be a linear combination of /i,/2 which is zero at v. Then is a harmonic 
function on Q which is zero on a biinfinite sequence of boundary points. We 
claim that this is impossible unless /a is identically zero. By the mean-value 
principle, from each zero of a (not identically zero) harmonic function on 
an infinite graph there are two paths to 00, one on which the function is 
increasing (and eventually positive) and another on which the function is 
decreasing (and eventually negative). Let a be a choice of increasing and 
eventually positive infinite path from a fixed lift v of v, and /3 be a choice of 
decreasing and eventually negative infinite path from the same point v. The 
union of the lifts of a starting at all lifts of is a nonnegative connected set 
of vertices containing all lifts of v, and containing a strictly positive path 
from —00 to 00; this contradicts the existence of f3. 

This completes the proof of the claim that /a must be identically zero, 
and thus the claim that dimWr = 1. 
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Fig. 2. A chain of loops. The roots o/detA are real and distinct. 

The roots of P vary continuously with the conductances. We can find 
a set of conductances for which all roots are real, see the next paragraph. 
As we vary the conductances, if two real roots ri,r2 collide and become 
a complex conjugate pair, the sum W^i © at the point of collision must 
be a two-dimensional subspace of the kernel of A, which we showed above 
cannot happen. Thus, the roots must remain real. From Theorem 1, P{zi) 
cannot have negative real roots so all roots are real and positive. 

We now find specific conductances for which all roots are real. Take a set 
of maximal cardinality of disjoint cycles, each winding around S, and com- 
plete it to a CRSF. Let Ci, . . . , be the cycles in order from innermost to 
outermost. Add k — 1 more edges to join these CRTs up into a connected 
set U (a "chain of loops," with bushes), so that Cj is connected to Cj+i with 
a path 7j. Put a large conductance R on every edge of U except for one edge 
of every Cj, and one edge of every path ji; these edges of U have weight 1. 
Put a small conductance e on the remaining edges of Q not in U. 

When £ is small and R is large, det A is to leading order a power of R 
times the determinant of the graph of Figure 2. This is because with high 
probability all edges with conductance R, none of the edges of conductance e, 
and some of the edges of conductance 1 will be present in a random CRSF. 
Contracting all edges of conductance R and removing those of conductance e, 
we are left with the graph of Figure 2 with all edges of conductance 1. For this 
graph, it is easy to verify that the roots of det A are real and distinct (see 
the example after Corollary 1). The actual Laplacian is a small perturbation 
of this so its determinant also has real distinct roots. □ 

Corollary 1. The number of cycles in a uniform random incompress- 
ible CRSF on an annulus is distributed as 1 plus a sum of k independent 
Bernoulli random variables, where k + 1 is the maximal number of nonin- 
tersecting incompressible cycles one can simultaneously draw on Q. 

Proof. Let w = 2 — z — 1/ z. Then Q{w) = P{z) has only real negative 
roots: 

k 

Q{w) = w\\{w + \i). 

i=l 

By Theorem 1, Q{w)/Q{l) is the probability generating function for the 
number of loops. It is also the probability generating function for a sum of 
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independent Bernoullis, with the ith Bernouhi being biased as {j^, i^^)- 
□ 

As an example, let us compute, for a rectangular cylinder, det A and 
the corresponding distribution of cycles in a uniform random incompress- 
ible CRSF. Let Qm be a line graph of length m and Hm,n = Gm x ^n- The 
eigenvalues of Ag^ are 2 + 2 cos ^ for k = 1,2, ... ,m [the corresponding 

eigenvectors are fk{x) = cos ^^^^^^^^^ for x = 0, 1, 2, . . . , m — 1] . By (3), we 
have 

detA = (2-z-z-^) Jj(^Ch„(A + 2) - z- 
= wYlw + Chn{X + 2)-2, 

A 

where w = 2 — z — ^ and Ch„ is defined by Ch„(Q; + ^) = a" + a""" (a variant 
of the Chebyshev polynomial). For nonnegative A and large n, Ch„(A + 2) 
is large unless A is near 0, so the relevant roots are A, = 2 + 2 cos for 
small j, j = 1, 2, .... We have 2 + A,- = 2 + ^ + 0(^)^ = a,- + 1/a,- where 
oij = 1 + ^ + Thus, in the limit m, n — )• oo with m/n — )• r the roots 

satisfy Ch„(Aj + 2) = 2 cosh 22 + o(l). 

The limit probability generating function for the number of cycles is then 

^/ \ /^/,\ /tf + 2cosh (vrj/r) — 2\ 

QM/Q(1) = . n ( 2cosh(vr,/r)- l ) ' 

For example, for a square annulus (m = n) the probability of a single cycle 
is approximately 95%. 

10.3.2. Torus example. We compute the Laplacian determinant for a flat 
line bundle on a n x n grid on a torus, that is, for the graph Z^/nZ^. In this 
case, the determinant is 

detA= n n 4-c-j-e-[ 

This product was evaluated in [1], yielding det A = C„P„(2;i, 2:2) where Cn 

is a constant tending to 00 and Pn tends to P, where 

(7) P{zi,Z2) = e-(^/2)(i2+fe2)^,^fc^_^^,+fc+,.fc_ 

We wish to rewrite this as 

Piz,,Z2) = E E C^kmi2 - Z{z\ - Z^=z:^^r, 

{j,k) primitive Tt^;! 
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where the first sum is over primitive vectors {j,k) in Z^, one per direction 
[i.e., only one of (j, k) and {—j, — k) appears]. Thus, Cjkm will be the relative 
probability that a CRSF has m cycles in homology class {j,k). To this end, 
pick a primitive vector (j, k) and consider the terms in (7) with monomials 
z^'' for md'L. Letting u = z-[z2 and q = e~^'^l'^^^^ ) the sum of these 
terms is 

which can be rewritten as a power series in 2 — u — 1/u as^ 

eez m=i \e=i ^ ^ / 

Thus, we find 

r , -V— r m + £ - l \ f Y+"^p-i^/^)0^+''^)e^ 
£=1 ^ ^ 
For example, the probability of a single cycle in homology class (1,0) is 

^-41%. 



10.3.3. Pants example. We consider here the case when S is a 3-holed 
sphere. This is the simplest case where the noncommutativity of the con- 
nection plays a role. We consider a flat SL2C-connection on {Q,T,). Let a,b 
and c be the monodromies around the three holes, which satisfy abc= 1. 
Since '7ri(S) is the free group on two generators, there are no other rela- 
tions. Note that a simple closed curve on S has only one of three possi- 
bly homotopy types: it must be homotopic to one of the three boundary 
curves. A finite lamination up to homotopy is therefore described by a triple 
{i,j,k) 7^ (0,0,0) of nonnegative integers, where there are i curves homo- 
topic to a, j homotopic to b, and k homotopic to c. Theorem 9 gives the 
following. 

Theorem 13. Let X = 2- Tr(a), F = 2 - Tr(6), Z = 2- Tr(c) = 2- 

Tr(a6). Then det = P{X,Y, Z) = Y^cksfs^'^^ '^^ere {i,j,k) is the 
number of cycles with homotopy type a,b,c, respectively. 



^We use the fact that 

r, i -I £ fm + £-l\, ,sm + l/r, 
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Note in particular that P is a polynomial with nonnegative coefficients. 
What polynomials occur? Is there an analog of Theorem 12 in this setting? 
We do not know. However, it is not difficult to show that the set of exponents 
{i,j,k) of the nonzero monomials of P{X,Y, Z) form a convex set in Z'^: 
let Pab be the length in terms of the number of vertices of the shortest 
path from face a (the face with monodromy a) to face b (the face with 
monodromy b). Similarly define pac and p^c- Then we can simultaneously 
embed i,j, k cycles of types a, b, c if and only if 



j + k< pbc- 

This is easily proved by induction on Pabi Pac, Pbc- So the Newton polyhedron 
of P is defined by i,j,k> 0, ii,j,k) / (0,0,0), and (8). 

11. Planar loop-erased random walk. In this section, we use the vector- 
bundle Laplacian to study the loop-erased walk on a planar graph. In Sec- 
tion 11.1, we compute the probability that the LERW between two boundary 
points passes left or right of a given face. In Section 11.2, we compute this 
probability for two faces. 

11.1. LERW around a hole. Let ^ be a finite planar graph. Let zi,Z2 
vertices on the outer boundary and / a bounded face. We compute the 
probability that the LERW from zi to Z2 goes left of /. This can be computed 
quite simply using duality: given a UST on a planar graph Q, the duals of 
the edges not in T form a UST of the planar dual graph. In the dual graph, 
our question is equivalent to asking whether the LERW from the face / first 
reaches the outer face along the edge between zi and Z2 or outside this edge. 
This is just the harmonic measure of the path from zi to Z2 as seen from /. 

Here we will compute this in another way, using line bundle technology. 
It illustrates a general method which will be useful in the next section. 

Take a line bundle with monodromy b around / and trivial around all 
other faces. This can be achieved by putting in a zipper of edges (dual of 
a simple path in the dual graph from / to the boundary) each with parallel 
transport b as in Figure 3. Suppose without loss of generality that the zipper 
ends between zi and Z2 as in the figure. 

Put in an extra edge from zi to Z2 with parallel transport 4>ziz2 = ^- Then 
by Theorem 1, 



(8) 



i+j< Pab, 
i + k<Pac, 



(9) 



detA=pi+p2[2-a--]+p3[2---- 




a I \ b a 
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b 



Fig. 3. Defining the connection for the LERW around a hole. 



Here p2 = P2{b) is the weighted sum of CRSFs which contain a cycle contain- 
ing edge a and going left of / (when oriented so that a is the edge before zi ) ; 
P3 =P3{b) is the weighted sum of CRSFs containing a cycle through a going 
"right" of /. 

We wish to compute p2{b = l),p3(6 = 1). These can be extracted from the 
coefficients of a and - in det A. 

a 

Letting zi,Z2 be the first two vertices, we have 

/ 1 -a 
A = A{b)+ -a-i 1 
V 0^ 

where A (6) is the Laplacian without the extra edge from zi to Z2 (and 
represents a vector or matrix of zeros). As a consequence, 

1 -a\ ( G{zi,zi) G{zi,Z2) 



det A = det A(6) det { I + \ ^ , , i i ^/ \ ^, \ 

V V~l/" 1 J \G{Z2,Zi) G(Z2,Z2) 

where G = A(6)"^ 

This last determinant is 1 + G{zi,zi) + G{z2,Z2) — aG{z2,zi) — MfliM^ 
Comparing coefficients with (9) yields 

b^Giz2,Zi)-G{z,,Z2) 



and 



P2 ■ 



P3 - 



52-1 

b{G{zi,Z2)-G{z2,Zi)) 
62-1 



detA(&) 



det A (6). 
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The desired quantity is 

P2 b^G{z2,zi) -G{zi,Z2) 



P2+P3 {b-l){bG{z2,Zi) + G{zi,Z2)) 

in the hmit 6 — )• 1. 

We need to compute G{v,v') when b is near 1. 

Lemma 2. When b = l + e, G{v,v') = ^^^^(^(l + eX{v,v') + 0{e'^)) 
where k is the number of spanning trees of Q and X{v,v') is the expected 
signed number of crossings of the zipper on a simple random walk from v 
to v'. 

Proof. For a matrix M, let denote (-l)*'4+«s times the deter- 
minant of M upon removing rows A and columns B. Here i^ is the sum 
of the indices of elements of A and likewise for is- (We can assume with- 
out loss of generality that all relevant vertices have even index so that we 
can ignore the signs.) Recall that for the standard Laplacian Aq, we have 
for any v,v' that (Aq)^ = n, the number of spanning trees of G. We also 
have (Ao)^^ = G{a,c) - G{a,d) - G{b,c) + G{b,d) where G is the Green's 
function.^ 

Letting Z = A{b) — Aq, a matrix supported on the zipper edges. We have 
G{v,v')detA{b) 

= A{b):, = (Ao + z):, = (Ao):, + zy {AoTz + • • • 

x,y 

(11) = + - l)(^o):;; + - l) (^o)X. + 0{e') 

uw 

= k(1 + £'^Go{v',u) + Go{v,w) 



-Go{v',w)-Go{v,u) + Oie^) 

Here the sum is over the zipper edges. For an edge uw the quantity Go(f u) + 
Go{v,w) — Go{v',w) — Go{v,u) is the transfer current, equivalently, the ex- 
pected signed number of crossings of uw of the simple random walk started 
at v' and stopped at v. □ 



®More generally, 

(Ao)^,V.v;>, = M«det(t + G(a„6,))i<. .<„ 

where [t]-P(t) denotes the coefficient of in P{t). This can be proved by induction on k 
and expanding the determinant on the left-hand side along row a^- 



SPANNING FORESTS AND THE VECTOR BUNDLE LAPLACIAN 



27 



By the lemma, when b = l + e the probabihty (10) becomes 

\ = IH 7, ^ 0{£) = 1 - X{zi,Z2) + 0(e) 

since X{zi,Z2) = -X{z2,zi). 

Corollary 2. The probability that the LERW from z\ to Z2 goes left 
of f is the total amount of current flowing left of f when 1 unit of current 
enters at z\ and exits at Z2. 

11.2. LERW around two holes. Let /i,/2 be two faces in a finite planar 
graph Q and zi,Z2 two points on the boundary. Consider the LERW from z\ 
to Z2. We compute the probability that it goes left of both faces. Although 
this can be computed in principle using Theorem 11 and the comments 
thereafter, in practice that method involves an integration of detA over 
the representation variety X; in the case of a large graph one might not 
have access to an explicit expression for detA. Here, we do the computation 
using only a bundle which is close to the trivial bundle; the determinant as 
a function on X can be expanded as a power series around the trivial bundle 
case and we only need the first two terms in the expansion, which we can 
write down explicitly in terms of the standard green's function. 

We use an SL2C-bundle. Choose a bundle with mondromy A around /i 
and B around /2, and trivial monodromy around the other faces. This can 
be achieved by putting trivial parallel transport on all edges except for two 
zippers, one from each of /i, /2 to the boundary, as in Figure 4. Let Za, Zb 
denote these two zippers. 

Add a new edge connecting zi to Z2 with parallel transport C. Suppose 
that A,B are close to the identity. Then a CRSF will have with high proba- 
bility only a loop containing C. The first-order correction to this consists in 




Fig. 4. Defining the connection for the LERW around two holes. 
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CRSFs having two loops, one of them containing C and the other of mon- 
odromy A^B or AB (or their inverses). The possible loops in a CRSF with 
one or two loops, one of which contains C, have one of the seven monodromy 
types (or their inverses): 

A, B, AB, C, CB^^A~\C{BA)''A^\BAy, C{ABy'B^\ABy, 

where n £ Z. 



Let 



A 



B 



C 



\ + xu y\fu 

z^/u 1 + wu 

1 + au h\fu 

Cy/u 1 + du 

e 1 
-1 



be the representations in SL2C, where u is small and x,y,z,w,a,b,c,d,e 
are variables. The equations det A = deti? = 1 yield x + w — yz = o{u) and 
a + d — bc = o{u). 

Computing the traces of the above seven classes, we have 

2 - Tr(^) = -{x + w)u + o{u), 

2 - Tr{B) = -{a + d)u + o{u), 

2 - Tt{AB) = -{a + d + w + X + cy + bz)u + o{u), 

(12) 2 - Tr(C) = 2 - e, 

2 - Tr(C(^S)"B-i(AS)-") = 2-e{l + du + o{u)), 

2 - TT{C{BA)''A~\BAy) = 2-e{l + wu + o{u)), 

2 - Tr{CB^^A^^) = 2-e{l + {d + w + bz)u). 

The relevant quantities for us are the last four and the products of any 
one the first three of these with the last four. 

We let A'^;; be the weight sum of configurations in which the LERW from zi 
to Z2 goes left of both faces, and having no other loops (i.e., CRSFs with one 
loop containing the extra edge, with the property that this loop surrounds 
both /i,/2). Similarly, define Nir,Nri,Nrr- Let Ni,N2,N3 be, respectively, 
the weighted sum of those configurations having one loop containing C and 
one loop surrounding, respectively, hole /i, hole /2, and both holes. 

If we extract the coefficient of e in — Qdet A, from (12) the result will be 
of the form 

Nil + Nir{l + du) + Nri{l + wu) + iV„.(l + {d + w + bz)u) 

- Ni{w + x)u - N2{a + d)u - N^^a + d + w + x + cy + bz)u + O(u^). 
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Recall that these variables are subject to the constraints x + w — yz = o{u) 
and a + d — bc = o(u). The coefficients of the A'^* are linearly independent 
given these constraints. Thus given an explicit expression for Qdet A to first 
order in u, it is a simple matter to extract these various coefficients A^*. The 
desired probabilities are 

Nil 

^" Na+Nir + Nrl + Nrr 

and so on. 

Let us show how to compute the coefficient of e in Qdet A, to first order 
in u. 

Indexing the vertices so that zi,Z2 are the first two vertices, we have 

/ / -CO' 
A{A,B,C) = A{A,B)+i-C-^ I 

V 0^ 

where A{A, B) is the Laplacian without the extra edge from zi to Z2 (and 
represents a vector or matrix of zeros). As a consequence, 

detA(A,S,C) 



detA(A,5)det / + 



/I -e -1^ 

110 

110 

V-1 -e 1 




ni (2) n( (2) (2) 

where G = A{A,B)~'^ . Note that G is self-dual, since it is the inverse of 
a self-dual matrix. Thus, we have 

Qdet A{A, B, C) = Qdet A(A, B) 

V ^ / VG(^1'^2) G{Z2,Z2) 

where G(zi,zi) and (js[z2,Z2) are diagonal matrices. The right-hand side is 
Qdet A ( A, times 

1 + - G{z^\zf^)+G{zf\z^^^) + G(4^.«) - eG(zf ). 

Thus, the coefficient of e in — Qdet A{A, B, C) is 



G{z[^\z^^^) Qdet A{A,B) 



to first order in u. 
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It is possible to give an exact expression for G{z\ )QdetA{A,B) as 
a sum over and Zb of products of standard Green's functions. Since A, B 
are close to the identity, we can write 

where Z is a perturbation supported on the zippers. As a consequence, both 
) , ) det AiA, B) = det A(A, B) 

and Qdet A(^,i?) itself have expansions in powers of u. The ratio of these 
is the desired quantity. 

(2') (2\ 

For the present purposes, we need the expansion of G{z\ , ) det A(^, B) 
to first order (which is the term of order ) and ^J^LetK{A^B) to first order 
(which is the term of order u). Unfortunately these computations, although 
elementary, can be quite long, and we have not been able to get an explicit 
answer for a main case of interest which is the upper half plane. 



12. Green's function and random walks. We give here a probabilistic 
interpretation of the Green's function of the line bundle Laplacian in the case 
where each edge has parallel transport either 1 or z. Let ^ be a graph, Ez 
a collection of directed edges, and consider the line bundle with parallel 
transport z on Ez and 1 elsewhere. Let G be the inverse of the line bundle 
Laplacian. 

Let v^v' be vertices of Q. Let us first ask: what is the distribution of the 
signed number of crossings of Ez on the SRW from v to v'l 

Let P be the transition matrix of the simple random walk on Q (ignoring 
the line bundle) with absorbing state v' . Thus, Pv2Vi is the probability of 
going to V2 from state vi (as usual the indices are reversed so that composi- 
tion respects the natural path order). Multiply the entry Py^v^ by z or I/2 if 
V1V2 is in Ez or —Ez-, respectively. With v' as the vertex with last index P 
has the form P = ( ^ ^ ) . Then the probability generating function of the 
number of crossings is 

00 

Y,{AN'').>. = {A{I-N)-\,^. 

k=0 



Proposition 3. The probability generating function of the signed num- 
ber of crossings of Ez on a simple random walk from v to v' is 
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Proof. Note that ((I - P)D)v^v^ = A^^^j for V2 / f', where D is the 
diagonal matrix of vertex degrees. As a consequence 



{A{I - iV)"i),,, = A,.,, {I - N)-i 

= 'YAy 



cof (J - A^)y 
'"det(/-iV) 



cof(AD"i):^r 



cof A;;^r'(deg(^;0deg(i;2)/detZ)) 



V'V2 



V2 



(det A)G(t;', t;Odeg(i;')/det D 



(detA)GK,.0?^^'^^"^^^"' 



G(z;',?;') 



□ 



The quantity 1 — ^z''^) jgg^' '^^^ ^ probabihstic interpretation. Choos- 
ing a gauge such that (p^i^ = 1 for all adjacent to v' , and using 

^ A„,„2G('U2,u') = Syy 

we have (setting v = v') 

degv'G{v\v')- G{v2,v') = l 

giving 

degf' ^ degu' ' 
Thus, we have the following proposition. 

Proposition 4. 1 — ^/-^ ^, is ^/le probability generating function 

of the signed number of crossings of Ez on a SRW started at v' and stopped 
on its first return to v' . 

Now let z = l + e with small e and 

(13) §^ = l + ^i^ + ^2S^ + ---, 

G{v' ,v') 
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so that Xi is the expected signed number of crossings of Ez on a SRW 
from V to v' . Similarly let 

(14) n( > \a 7 = l + ^2e' + l3e' + ---, 

G{v',v') deg v' 

so that 212 is the expectation of the square of the signed number of crossings 
of Ez of the SRW from v' until its first return to v' (note that Yi is zero 
and = — 12 by symmetry). 

Proposition 5. The determinant of the line bundle Laplacian A above 
satisfies 

(15) det A = Kdegv'{-Y2e'^ + l2e^ + 0{e^)), 

where 2Y2 is the expected square of the signed number of crossings of Ez of 
a SRW from v' until its first return to v' . In particular Y2degv' does not 
depend on v' . 

Proof. From Lemma 2, we have 

G{v, v') det A = k(1 + Xie + 0{e^)). 
However by (13) and (14), we have 

n( '^A,A 1 l + Xi£ + ^£2+0(£3) 

G[v,v jdet A = — -—5 — — det A, 

so we must have det A = Kdegf'(— 12^^ + ¥28^ + 0(e^)). □ 

We can get a full expansion of det A as a power series in z around z = 1 
as follows. We have 

^logdetA= Gl-^^Gl. 

uwdEz 

This follows from differentiating the entries in A. Here using Propositions 4 
and 3, the quantities on the right-hand side have explicit probabilistic in- 
terpretations. Integrating both sides from z = 1 to z = zq, gives the desired 
expansion. 

13. Monotone lattice paths. Let Tm,n be the mx n grid on a torus [the 
nearest neighbor grid {0, 1, ... , m} x {0, 1, . . . , n} with opposite sides iden- 
tified: {m,j) ~ (0,j) and {i,n) ~ («,0)]. Choose for each vertex a north- or 
east-going edge, independently and with probability p= 1/2. The resulting 
configuration of edges makes a directed CRSF, and we wish to determine the 
distribution of the number and homology type of the cycles. This problem 
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was studied in [7] who showed among other things that when m/n \s close 
to a rational with small denominator the number of cycles is (in the limit 
771,71—7-00) tending to a Gaussian with expectation on the order of ^/n. 

Here we show how to compute explicitly for each 777, 77 the probability gen- 
erating function of the total homology class of the cycles (which determines 
both the number of cycles and their direction). 

We make Tm,n into a directed graph, directing all edges northward or 
eastward. Let $ be a flat line bundle on Vm,n with monodromy z and w 
for loops in homology class (1,0) and (0, 1), respectively. <1> can be obtained 
by putting parallel transports 1 on all edges except edges (777 — 1,7)(0,7) 
which get parallel transport z and edges (7,77 — 1)(7,0) which get parallel 
transport w. 

By Theorem 6, the determinant of the line bundle Laplacian is 

(16) F^,n(z,7x;) = detA= Cj^til-z^w")', 

j,k>0 

where (j, k) = {p£, qi) and p, q are coprime, and Cj^k is the number of CRSFs 
with total homology class {j,k), that is, with i cycles of homology {p,q). 

This determinant can be explicitly computed using a standard Fourier 
diagonalization of A. 

Proposition 6. 

(17) Fm,n{z,w)=Yl YI2-U-V. 

u™- =z v"=w 

Proof. If u"^ = z and = w, then the function f{x,y) = u^v^ is an 
eigenvector of A with eigenvalue 2 — u — v. These eigenvectors are indepen- 
dent and span C^"" . □ 

From (17), we can extract the coefficients Cm,n- Let us consider for sim- 
plicity the case m = n. In [7], it was shown that with probability tending 
to 1 as 77 — )• 00, all cycles will have homology class (1, 1). Thus, from (16), 

(18) Fn,n{z,w)=J2CjA^-ZWy + ^ Cj^k{l - Z^w'^Y , 

j>0 j^k>0 

where the second sum is negligible, in the sense that the sum of its coefficients 
is o^'^Cjj). We can thus ignore the second sum, let z = 1 and expand 
Fn^ni^iw) around w = 1. We have 

Fn,n{z,w)= Yl 2-77-7;= (2-7i)"-77;. 

u"=z,v"=w u"=z 
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Letting Y = 1 — w and w) = Hn{l, Y), we have 

n-l 

Hn{i, y) = n (2 - e^-^/")" - (1 - y) 

i=o 

and 

i7„(l,l) ~ (2 - e2Tii/n)n • 

The product (2 — g^'^*-?/")" is large in absolute value unless \j\ is small, that 
is, when \j\ is 0(n^/^), so the dependence on Y comes from terms with \j\ 
small. That is, the product can be written as 

n (•■•) n (■•■) 

|j|<nl/2+e |j|>„l/2+e 

and the second product is 1 + o(e~"''). When \9\ is small, 

(2 - e*^)" = {l-ie + + •••)'' = e-*"^+'^^'+-. 
Plugging in ^ = 27rj'/n yields 

j=-oo 

where we have extended the range of j without loss of precision. Thus up 
to exponentially small errors, 

^n{l,Y) = y _ g-{27ri)2/n j_ yg-(27ri)2/n) 

is the probability generating function for the number -/V(x of (1, l)-cycles. 
We see that the number of (1,1) cycles is a sum of independent Bernoulli 
random variables with the jth having bias pj = e"^^'^-'^ Z"^. The expectation 
of A^(i,i) is then (using t = 27rj/-^/n) 




The variance is 



E\/n f°° 4-2, .2, , Fn f 1 \ 

Since the variance tends to oo, the central limit theorem implies that distri- 
bution tends to a Gaussian as n — )• oo. 

A similar computation holds when m = pk,n = qk and /c — )■ oo. 
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14. Open questions. 

1. Consider a graph embedded on a Riemann surface, sucli tliat edges are 
geodesic segments. Now use a naturally associated connection, for ex- 
ample, the Levi-Civita connection on the tangent bundle. What are the 
probabilistic consequences of choosing such a connection? 

2. Is there any combinatorial meaning to the vector bundle Laplacian for 
higher-rank bundles? 

3. What information does the line bundle Laplacian give about the three- 
dimensional monotone lattice path problem, generalizing the results of 
Section 13? 
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